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PREFACE 

This problem was born in 1922 when A.Friedmann wrote 
his famous cosmological solution for the homogeneous 
isotropic Universe. However, during the next 35 years 
researches devoted their attention mainly to the physi- 
cal processes after the Big Bang and there was no se- 
rious attempts to put under a rigorous analysis the phe- 
nomenon of the cosmological singularity as such. The 
person who inspired the beginning of such analysis was 
L.D. Landau. In the late 1950s he formulated the crucial 
question whether the cosmological singularity is a gen- 
eral phenomenon of General Relativity or it appears only 
in particular solutions under the special symmetry con- 
ditions. The large amount of work have been done in 
Landau school before an answer emerges in 1969 in the 
form of the so-called "BKL conjecture" (the present-day 
terminology). The basic reviews 1 covering also the con- 
temporary development are [ 1 ] - [ 1 4] . The BKL conjec- 
ture has its foundation in a collection of results and, first 
of all, it asserts that the general solution containing the 
cosmological singularity exists. This fundamental ques- 
tion of existence of such solution was the principal goal 



1 Some of these articles should not be considered merely as reviews 
since they contain also the original results. This is relevant especially 
for the papers [3], [8] and [13]. For example, the first construction of the 
general non oscillatory solution with power asymptotic near singularity 
for the case of perfect liquid with stiff matter equation of state and 
investigation of the general oscillatory regime in presence of the SU(2) 
Yang-Mills field can be found only in [3]. The Iwasawa decomposition 
for the metric tensor was first introduced namely in [8] . Many original 
details of the dynamical system approach to cosmological evolution 
can be found in [13]. The papers [13] and [14] we mentioned for 
completeness in order to remind to a reader on the existence of an 
alternative approach to the analysis of the character of cosmological 
singularity based on the representation of the Einstein equations in 
the form of dynamical system and the serch for the description of its 
attractor in vicinity to the singularity. This is powerful method which 
can be considered as dual to the "cosmological billird approch" to 
which the present talk is mainly dedicated. The restricted time for the 
talk prevented to include a discussion also on this important aspect of 
the theory. 



of our work, however, we succeeded also in describing 
the analytical structure of gravitational and matter fields 
in asymptotic vicinity to the singularity and we showed 
that in most general physical settings such solution has 
complicated oscillatory behaviour of chaotic character. 

In order to avoid misunderstandings let's stress that 
under cosmological singularity we mean the singularity 
in time, when singular manifold is space-like, and when 
the curvature invariants together with invariant character- 
istics of matter fields (like energy density) diverge on this 
manifold. An intuitive feeling that there are no reasons 
to doubt in existence of the general solution with cosmo- 
logical singularity we have already in 1964 but another 
five years passed before the concrete structure have been 
discovered. In 1965 appeared the important theorem of 
Roger Penrose [15], saying that under some conditions 
the appearance of incomplete geodesies in space-time is 
unavoidable. This is also singularity but of different type 
since, in general, incompleteness does not means that in- 
variants diverge. In addition the theorem can say noth- 
ing about the analytical structure of the fields near the 
points where geodesies terminate. Then Penrose's result 
was not of a direct help for us, nevertheless it stimu- 
lated our search. Today it is reasonable to consider that 
the BKL conjecture and Penrose theorem represent two 
sides of the phenomenon but the links are still far to be 
understandable. This is because BKL approach deal with 
asymptotic in the vicinity to the singularity and Penrose 
theorem has to do with global space-time. 

It is worth to stress that some misleading statements 
are to be found in the literature in relation to the afore- 
mentioned results. First of all, from the BKL theory as 
well as from Penrose theorem not yet follows that cos- 
mological singularity is inevitable in General Relativity. 
BKL showed that the general solution containing such 
singularity exists but general in the sense that initial data 
under which the cosmological singularity is bound to ap- 
pear represent a set of nonzero measure in the space of 
all possible data. However, we don't know "how big" this 
measure is and we have no proof that this set can cover 



the totality of initial data. In the non-linear system can be 
many general solutions (that is, each containing maximal 
number of arbitrary functional parameters) of different 
types including also a general solution without singular- 
ity. Moreover there is the proof [16] of the global sta- 
bility of Minkowski spacetime which means that at least 
in some small (but finite) neighbourhood around it exists 
a general solution without any singularity at any time. 
The same is true in relation to the all versions of Penrose 
theorem: for these theorems to be applicable the nontriv- 
ial initial/boundary conditions are strictly essential to be 
satisfied, but an infinity of solutions can exists which do 
not meet such conditions. The thorough investigation of 
applicability of Penrose theorem as well as all its subse- 
quent variations the reader can find in [17, 18]. 

The second delusion is that the general solution with 
singularity can be equally applied both to the singular- 
ity in future (Big Crunch) and to the singularity in past 
(Big Bang) ignoring the fact that these two situations 
are quite different physically. To describe what are go- 
ing near cosmological Big Crunch (as well as near the 
final stage of gravitational collapse of an isolated object 
in its co-moving system) one really need the general so- 
lution since in the course of evolution inescapably will 
arise the arbitrary perturbations and these will reorganize 
any regime into the general one. The Big Bang is not the 
same phenomenon. We don't know initial conditions at 
the instant of singularity in principle and there are no rea- 
sons to expect that they should be taken to be arbitrary. 
For example, we can not ruled out the possibility that the 
Universe started exactly with the aid of the Friedmann 
solution and it may be true that this does not means any 
fine tuning from the point of view of the still unknown 
physics near such exotic state. Of course, the arbitrary 
perturbations familiar from the present day physics will 
appear after Big Bang but this is another story. The con- 
clusion is that if somebody found the general cosmolog- 
ical solution this not yet means that he knows how Uni- 
verse really started, however he has grounds to think that 
he knows at least something about its end. 

Sometimes one can find in literature the statement that 
in the BKL approach only the time derivatives are impor- 
tant near singularity and because of this the asymptotic 
form of Einstein equations became the ordinary differen- 
tial equations with respect to time. Such statement is a lit- 
tle bit misleading since space-like gradients play the cru- 
cial role in appearing the oscillatory regime. One of the 
main technical advantage of the BKL approach consists 
in identification among the huge number of the space 
gradients those terms which are of the same importance 
as time derivatives. In the vicinity to the singularity these 
terms in no way are negligible, they act during the whole 
course of evolution (although from time to time and dur- 
ing comparatively short periods) and namely due to them 
oscillations arise. The subtle point here is that asymp- 



totically these terms can be represented as products of 
effective scale coefficients, governing the time evolu- 
tion of the metric, and some factors containing space- 
like derivatives. This nontrivial separation springing up 
in the vicinity to the singular point produce gravitational 
equation of motion which effectively are the ordinary dif- 
ferential equations in time because all factors containing 
space-like derivatives enter these equations solely as ex- 
ternal parameters, though dynamically influent parame- 
ters. Owing to these ordinary differential equations the 
asymptotic evolution can be described as motion of a 
particle in some external potential. The aforementioned 
dominating space gradients create the reflecting potential 
walls responsible for the oscillatory regime. For the case 
of homogeneous cosmological model of the Bianchi IX 
type such potential have been described by Misner [19]. 
The literal assertion that "only the time derivatives are 
important near singularity" is correct just for those cases 
when the general solution is of non oscillatory character 
and has simple power asymptotic near singularity as, for 
instance, for the cases of perfect liquid with stiff-matter 
equation of state [3, 20, 21], pure gravity in space-time 
of dimension more than ten [22], or some other classes 
of "subcritical" field models [23]. 

BASIC STRUCTURE OF 
COSMOLOGICAL SINGULARITY 

The character of the general cosmological solution in 
the vicinity to the singularity can most conveniently be 
described in the synchronous reference system, where 
the interval is of the form 

-ds 2 = -dt 2 +g afj dx a dx P (1) 

We use a system of units where the Einstein grav- 
itational constant and the velocity of light are equal 
to unity. The Greek indices refer to three-dimensional 
space and assume the values 1,2,3. Latin indices i,k 
will refer to four-dimensional space-time and will take 
the values 0,1,2,3. The coordinates are designated as 
(x° ,x l ,x 2 ,x 3 ) = {t,x,y,z). 

The Einstein equations in this reference system take 
the form 

R° = l -k+ l -K$4=T° - l -T, (2) 
R P a = ^(V84)+P^ = T P a-^T , (4) 



where the dot signifies differentiation with respect to 
time t and 

K a p = g a p , 8 = detg aj3 . (5) 

The tensorial operations on the Greek indices, as well 
as covariant differentiation in this system are performed 
with respect to the three-dimensional metric g a p. The 
quantity K is a three-dimensional contraction: 

K=K« = Qng). (6) 

Pa is a three-dimensional Ricci tensor, expressed in 
terms of g a p in the same way as PJf is expressed in terms 

of gfo. The quantities T q, T {) a and T a are components of 
the energy-momentum tensor T * four-dimensional con- 
traction of which is designated by T . 

T =t£ = T q + T ® . (7) 

It turn out that the general cosmological solution of 
Eqs. (2)-(4) in the asymptotic vicinity of a singularity 
with respect to time is of an oscillatory nature and may 
be described by an infinite alternation of the so-called 
Kasner epochs. The notions of a Kasner epoch and of the 
succession of two of these epochs are the key elements 
in the dynamics of the oscillatory regime. It is most 
convenient to study their properties in the example of 
empty space, when T * = and then take into account 
all the changes that may be observed in the presence of 
matter. This procedure is reasonable since, in general, the 
influence of matter upon the solution in the vicinity of the 
singularity appears to be either negligible or can be put 
under the control. 

So let's assume that the tensor T \ in Eqs.(2)-(4) 
equals zero. A Kasner epoch is a time interval during 
which in Eq. (4) the three-dimensional Ricci tensor P„ 
may be neglected in comparison with the terms involving 
time differentiation. Then from (2) and (4) we obtain the 
following equations in this approximation: 

(Vg4) = o, = o. (8) 

Here and elsewhere we shall assert that the singularity 
corresponds to the instant t = and we shall follow the 
evolution of the solution towards the singularity, i.e., the 
variation of time as it decreases from certain values t > 
down to t = 0. Eq. (3) in the general case is of no interest 
for the dynamics of the solution, since its role is reduced 
to the establishment of certain additional relations on 
arbitrary three-dimensional functions resulting from the 
integration of the Eqs. (2), (4) (that is of supplementary 
conditions for initial data). 

The general solution of Eqs. (8) may be written down 
in the form 



Sa/»=W£/f, T lAB =diag(t 2 P\t 2 Pi,t 2 Pi) (9) 

where by the big Latin letters A,B,C we designate the 
three-dimensional frame indices (they take the values 
1,2,3). The exponents pi,P2,P3 and vectors 1^ are arbi- 
trary functions of the three-dimensional coordinates x a . 
We call the directions along 1^ as Kasner axis, the triad 
la represents the common eigenvectors both for metric 
g a p and second form K a p. The exponents p\,P2,P3 sat- 
isfy two relations: 

Pl+P2+P3 = l, p\+P2 + p\ = \- (10) 
It ensues from these relations that one of the expo- 
nents pa is always negative while the two others are pos- 
itive. Consequently, the space expands in one direction 
and contracts in two others. Then the value of any three- 
dimensional volume element decreases since, according 
to (9)-(10), the determinant of g a p decreases proportion- 
ally to f 2 . 

Of course, the solution (9)-(10) sooner or later will 
cease to be valid because the three-dimensional Ricci 
tensor P„ contain some terms which are growing with 
decreasing of time faster than the terms with time deriva- 
tives and our assumption that P„ can be neglected will 
become wrong. It is possible to identify these "danger- 
ous" terms in P„ and include them into the new first 
approximation to the Einstein equations, instead of (8). 
The remarkable fact is that the asymptotic solution of this 
new approximate system can be described in full details 
and this description is valid and stable up to the singu- 
larity. The result is that the evolution to the singularity 
can be represented by a never-ending sequence of Kas- 
ner epochs and the singularity t = is the point of its 
condensation. The durations of epochs tend to zero and 
transitions between them are very short comparatively to 
its durations. The determinant of the metric tensor g a p 
tends to zero. On each Kasner epoch the solution take 
the form (9)-(10) but each time with new functional pa- 
rameters pa and l^. On each epoch the exponents p A sat- 
isfy the same relations (10), that is the space expands in 
one direction and contracts in two others, however, from 
epoch to epoch these directions are different, i.e. on each 
new epoch the Kasner axis rotate relatively to their ar- 
rangement at the preceding one. 

The effect of rotation of Kasner axis make its use in- 
convenient for an analytical description of the asymp- 
totic oscillatory regime because this rotation never stops. 
However, it turn out that another axis exist (they are not 
eigenvectors for the second form K a p) , rotation of which 
are coming to stop in the limit t — > and projection of 
the metric tensor into such "asymptotically frozen" (ter- 
minology of the authors of Ref. [8]) triad still is a di- 
agonal matrix. The components of this matrix have no 



limit since their behaviour again can be described by 
the never-ending oscillations of a particle against some 
potential walls. This is an efficient way to reduce the 
description of asymptotic evolution of six components 
of the metric tensor to the three oscillating degrees of 
freedom. For the homogeneous model of Bianchi type 
IX this approach was developed in [24, 25] where the 
three-dimensional interval has been represented in the 
form g ap dx a dxl i = {RTR) AB (l^dx a )(lfdxP) with the 

standard Bianchi IX differential forms l A dx a (where l A 
depends only on x a in that special way that d v l^ — 
— Cab fa v with on ty non- vanishing structural con- 
stants C23 = Cfj = C 3 2 = !)■ The diagonal matrix T and 
three-dimensional orthogonal matrix R depend only on 
time (tilde means transposition). Remarkably, the gravi- 
tational equations for this model shows that near singu- 
larity all three Euler angles of matrix R tends to some 
arbitrary limiting constants and three components of T 
oscillate between the walls of potential of some special 
structure. 

We never tried to generalize this approach (namely 
with orthogonal matrix R) to the inhomogeneous mod- 
els but the recent development of the theory showed 
that even in most general inhomogeneous cases (includ- 
ing multidimensional spacetime filled by different kind 
of matter) there is analogous representation of the met- 
ric tensor leading to the same asymptotic freezing phe- 
nomenon of "non-diagonal" degrees of freedom and re- 
ducing the full dynamics to the few "diagonal" oscillat- 
ing scale factors. This is so-called Iwasawa decomposi- 
tion first used in [8] and thoroughly investigated in [26]. 
The difference is that in general inhomogeneous case in- 
stead of orthogonal matrix R it is more convenient to 
use an upper triangular matrix N (with components N„ 
where upper index A numerates the rows and lower index 
a corresponds to columns) 



space) which follow from the Lagrangian: 



N = 




(11) 



and to write three-dimensional interval in the form 
g aj} dx a dxP = (NTN) a pdx a dxP). The diagonal matrix 
r as well as matrix N are functions of all four coordi- 
nates but near the singularity matrix N tends to some 
time-independent limit and components of F oscillate be- 
tween the walls of some potential. This asymptotic os- 
cillatory regime has the well defined Lagrangian. If one 
writes matrix T as T = diag(e^ 2 P i ,e~ 2 P 2 ,e~ 2 P 3 ) then 
the asymptotic equations of motion for the scale coef- 
ficients j3 A became the ordinary differential equations in 
time (separately for each point x a of three-dimensional 



T ~ dp A dp B 
L - (j AB — - dT 



v(P A ) 



de^P 1 + C 2 e-W 2 -^ +C,e-W i -f il ) . 



(12) 

Here we use the new time variable T instead of original 
synchronous time t. In asymptotic vicinity to the singu- 
larity the link is dt = — yj&etg a pd% where differentials 
should be understood only with respect to time, consider- 
ing the coordinates x a in detg a p formally as fixed quan- 
tities. Since detg a p tends to zero approximately like t 2 it 
follows that singular limit t — > corresponds to T — > °°. 
The metric Gab of three-dimensional space of scale co- 
efficients j3 A are defined by the relation GABdfi A df5 B = 
L(</j3 A ) 2 - (Zdfi A ) 2 - This is flat Lorenzian metric with 
signature (-,+,+) which can be seen from transformation 
j3' =j3 1 +/3 2 + j3 3 ,/3 2 = j3 1 -j3 2 + j3 3 , j3 3 = -j3 3 af- 
ter which one get G AB dl5 A dl5 B = -2(dj3 l ) 2 +2(dj3 2 ) 2 + 
2(<fj3 3 ) 2 . All coefficients Ca(x") are time-independent 
and positive; with respect to the dynamics they play a 
role of external fixed parameters. Apart from the three 
differential equations of second order for j3 A which fol- 
low from the Lagrangian (12) there is well known addi- 
tional constraint 



dH A dl3 B A 

Gab 7777 + v(/3 ) = ' 



(13) 



which represents the (q) component of the Einstein 
equations. In particular case of homogeneous model of 
Bianchi type IX equations (12) and (13) gives exactly 
the same system which was described in [24, 25], in 
spite of the fact that in these last papers the asymptot- 
ical freezing of "non-diagonal" metric components has 
been obtained using an orthogonal matrix R instead of 
Iwasawa's one N. Analysis of the eqs. (12)-(13) shows 
that in the limit T — > °° the exponents j3 A (r) are positive 
and all tend to infinity in such a way that the differences 
/3 2 — j3 1 and j3 3 — j3 2 also are positive and tend to infin- 
ity, that is each term in the potential V(j3 A ) tends to zero. 
Then from (13) follows that each trajectory )3 a (t) be- 
comes "time-like" with respect to the metric Gab, i- e - 

near singularity we have Gab-^-^ < 0, though this 
is so only in the extreme vicinity to the potential walls 
j3 1 = 0, j3 2 - j3 1 = and j3 3 - j3 2 = 0. Between the walls 
where j3 1 > 0, j3 2 - j3 1 > 0, j3 3 - j3 2 > the potential is 
exponentially small and trajectories become "light-like", 

i.e. Gab^-^jT = 0. These periods of "light-like" mo- 
tion between the walls corresponds exactly to the Kas- 
ner epochs (9)-(10) (with an appropriate identification 
of Kasner axis during each period). It is easy to see that 
the walls itself are "time-like" what means that collisions 
of a "particle moving in a light-like directions" against 
the walls are inescapable and interminable. 



One of the crucial points discovered in [8] is that in the 
limit T — > °° the walls become infinitely sharp and of infi- 
nite height which simplify further the asymptotic picture 
and make transparent the reasons of chaoticity of such 
oscillatory dynamics. Because GabP A P B = — 2(/3 1 j3 2 + 
j3 1 j8 3 + J3 2 j3 3 ) and near singularity all j3 A are positive we 
have G A bP A P B < 0. Then by the transformation fi A = 
pY^ one can introduce instead of j3 A the "radial" coordi- 
nate p > (p — > °° when T — > °°) and "angular" coordi- 
nates y 4 subjected to the restriction Gab^J 8 = — 1- The 
last condition pick out in y-space the two-dimensional 
Lobachevsky surface of constant negative curvature and 
each trajectory J3 (t) has the radially projected trace on 
this surface. The free Kasner flights in three-dimensional 
/3-space between the walls are projected into geodesies 
of this two-dimensional surface. The walls are projected 
into three curves forming a triangle on the Lobachevsky 
surface and reflections against these curves are geometric 
(specular). If we introduce the new evolution parameter 
T by the relation dx = p 2 dT then the new Lagrangian 
(with respect to the "time" T) will be: 

LT = -(^) 2 + G A B^-p 2 V(pf), 
G AB ff = -l- 

(14) 

In the limit p = °° the new potential p 2 V(py A ) is 
exactly zero in the region between the walls where y 1 > 
0, y 2 — y 1 > 0, y 3 — y 2 > and becomes infinitely large 
at the points y 1 = 0, y 2 - y 1 = 0, y 3 — y 2 = where 
the walls are located and behind them where quantities 
y\ y 2 — y',y 3 — y 2 are negative. This means that near 
singularity potential p 2 V depends only on y- variables 
and p can be considered as cyclic degree of freedom. In 
this way the asymptotic oscillatory regime can be viewed 
as the eternal motion of a particle inside a triangular 
bounded by the three stationary walls of infinite height 
in two-dimensional space of constant negative curvature. 
The important fact is that the area occupied by this 
triangle is finite. It is well known (see references in 
[8], section 5.2.2) that the geodesic motion under the 
conditions described is chaotic. 

It is worth to mention that in case of homogeneous 
Bianchi IX model the fact that its dynamics is equivalent 
to a billiard on the Lobachevsky plane was established in 
[27]. 

The numerical calculations confirming the admissi- 
bility of the BKL conjecture can be found in [6] and 
[28,29]. 



THE INFLUENCE OF MATTER 

In papers [3, 20, 30] we studied the problem of the in- 
fluence of various kinds of matter upon the behaviour 



of the general solution of the gravitational equations in 
the neighbourhood of a singular point. It is clear that, de- 
pending on the form of the energy-momentum tensor, we 
may meet three different possibilities: (i) the oscillatory 
regime remains as it is in vacuum, i.e. the influence of 
matter may be ignored in the first approximation; (ii) the 
presence of matter makes the existence of Kasner epochs 
near a singular point impossible; (iii) Kasner epochs ex- 
ist as before, but matter strongly affects the process of 
their formation and alternation. Actually, all these possi- 
bilities may be realized. 

There is a case in which the oscillatory regime ob- 
served as a singular point is approached remains the 
same, in the first approximation, as in vacuum. This case 
is realized in a space filled with a perfect liquid with the 
equation of state p = ke for < k < I. No additional re- 
flecting walls arise from the energy-momentum tensor in 
this case. 

If k — 1 we have the "stiff matter" equation of state 
p = e. This is the second of the above-mentioned possi- 
bilities when neither Kasner epoch nor oscillatory regime 
can exist in the vicinity of a singular point. This case 
has been investigated in [3, 20] where it has been shown 
that the influence of the "stiff matter" (equivalent to the 
massless scalar field) results in the violation of the Kas- 
ner relations (10) for the asymptotic exponents. Instead 
we have 

Pi+P2 + P3 = h p\+pI + pI = 1-P% (15) 

where p 2 ^ is an arbitrary three-dimensional function 
(with the restriction p 2 ^ < 1) to which the energy den- 
sity e of the matter is proportional (in that particular 
case when the stiff-matter source is realized as a mass- 
less scalar field <p its asymptotic is (p = p^lnt and this is 
the formal reason why we use the index <p for the addi- 
tional exponent p,p). 

Thanks to (15), in contrast to the Kasner relations 
(10), it is possible for all three exponents p A to be pos- 
itive. In [20] it has been shown that, even if the con- 
traction of space starts with the quasi-Kasner epoch (15) 
during which one of the exponents p A is negative, the 
asymptotic behaviour (9) with all positive exponents is 
inevitably established after a finite number of oscillations 
and remains unchanged up to the singular point. Thus, 
for the equation of state p = e the collapse in the gen- 
eral solution is described by monotonic (but anisotropic) 
contraction of space along all directions. The asymp- 
totic of the general solution near cosmological singular- 
ity for this case we constructed explicitly in [3], see also 
[20, 21, 23]. The disappearance of oscillations for the 
case of a massless scalar field should be consider as an 
isolated phenomenon which is unstable with respect to 
inclusion into the right hand side of the Einstein equa- 
tions another kind of fields. For instance, in the same pa- 



per [20] we showed that if to the scalar field we add a 
vector one then the endless oscillations reappear. 

The cosmological evolution in the presence of an elec- 
tromagnetic field may serve as an example of the third 
possibility. In this case the oscillatory regime in the pres- 
ence of matter is, as usual, described by the alternation of 
Kasner epochs, but in this process the energy-momentum 
tensor plays a role as important as the three-dimensional 
curvature tensor. This problem has been treated by us 
in [30], where it has been shown that in addition to the 
vacuum reflecting walls also the new walls arise caused 
by the energy-momentum tensor of the electromagnetic 
field. The electromagnetic type of alternation of epochs, 
however, qualitatively takes place according to the same 
laws as in vacuum. 

In paper [3] we have also studied the problem of 
the influence of the Yang-Mills fields on the charac- 
ter of the cosmological singularity. For definiteness, we 
have restricted ourselves to fields corresponding to the 
gauge group SU(2). The study was performed in the syn- 
chronous reference system in the gauge when the time 
components of all three vector fields are equal to zero. 
It was shown that, in the neighbourhood of a cosmolog- 
ical singularity, the behaviour of the Yang-Mills fields is 
largely similar to the behaviour of the electromagnetic 
field: as before, there appears an oscillatory regime de- 
scribed by the alternation of Kasner epochs, which is 
caused either by the three-dimensional curvature or by 
the energy-momentum tensor. If, in the process of al- 
ternation of epochs, the energy-momentum tensor of the 
gauge fields is dominating, the qualitative behaviour of 
the solution during the epochs and in the transition re- 
gion between them is like the behavior in the case of free 
Yang-Mills fields (with the Abelian group). This does 
not mean that non-linear terms of the interaction may be 
neglected completely, but the latter introduce only mi- 
nor, unimportant quantitative changes into the picture we 
would observe in the case of non-interacting fields. The 
reason for this lies in the absence of time derivatives of 
the gauge field strengths in those terms of the equations 
of motion which describe the interaction. 



MULTIDIMENSIONAL SPACETIME 
AND SUPERGRAVITY 

The story resembling the aforementioned effect of dis- 
sappearance (for scalar field) and reconstruction (after 
adding a vector field) of oscillations occurred later in 
more general and quite different circumstances. In 1985 
appeared very interesting and unexpected result [22] 
that oscillatory regime near cosmological singularity in 
multidimensional spacetime (for pure gravity) holds for 
spacetime dimension D up to D = 10 but for dimension 



D > 1 1 the asymptotic of the general solution follow the 
smooth multidimensional Kasner power law. Up to now 
we have no idea why this separating border coincides 
with dimension so significant for superstring theories, 
most likely it is just an accident. However, the important 
point is that if we will add to the vacuum multidimen- 
sional gravity the fields of /?-forms the presence of which 
is dictated by the low energy limit of superstring mod- 
els, the oscillatory regime will reappear. This fact was 
established in [31, 32] and subsequently has been devel- 
oped by TDamour, M.Henneaux, H.Nicolai, B Julia and 
their collaborates into the new interesting and promis- 
ing branch of superstring theories. In articles [31, 32] it 
was demonstrated that bosonic sectors of supergravities 
emerging in the low energy limit from all types of super- 
string models have oscillatory cosmological singularity 
of the BKL character. Let consider the action of the fol- 
lowing general form: 



2^ 



(p+l)\ 



R - d'(pdj(p- 

e X p <Pp(P+l) p(p+l)h-i p +i 
h—'p+i 



(16) 



where designates the the field strengths generated 

by the /?-forms A p , i.e. F^ 1 ^ = atitisym(d il A i2 ... ip+l ). 

The real parameters X p are coupling constants corre- 
sponding to the interaction between the dilaton and 
/?-forms. The tensorial operations in (16) are earring 
out with respect to D-dimensional metric g% and 
g = |detgHt|. Now the small Latin indices refer to D- 
dimensional space-time and Greek indices (as well as 
big Latin frame indices A,B and C) correspond to d- 
dimensional space where d = D — 1 . Also in this theory 
the Kasner-like epochs exist which are of the form: 



g ik dytdx k = -dt 2 + ri AB (t,x a )lf l (x a )^{x a )dx^dx v , 
ri AB = diag[t 2 P^ a \t 2 P^ a \...,t 2 P^ x % 



<p = p v (x a )lnt + <po(x a ). 



(17) 
(18) 



However, in the presence of the dilaton the exponents p A 
instead of the Kasner law satisfy the relations analogous 
to (15): 

X>a = 1, Zp 2 a = 1-pI (19) 



A=l 



The approximate solution (17)-(19) follows from the 
D-dimensional Einstein equations by neglecting the 
energy-momentum tensor of /5-forms, J-dimensional 
curvature tensor and spatial derivatives of (p. Now 
one has to do the work analogous to that one for 4- 
dimensional gravity: it is necessary to identify in all ne- 
glected parts of the equations those "dangerous" terms 
which will destroy the solution (17)-(19) in the limit 



t — > 0. Then one should construct the new first approx- 
imation to the equations taking into account also these 
"dangerous" terms and try to find asymptotic solution 
for this new system. This is the same method which have 
been used in case of the 4-dimensional gravity with elec- 
tromagnetic field and it works well also here. Using the 
Iwasawa decomposition for <i-dimensional frame metric 
r\ AB = (NTN) AB where T = diag (e~ 2 ^ ' , e' 2 ^ , . . . , e~ 2 ^ ) 
it can be shown [8] that near singularity again the phe- 
nomenon of freezing of "non-diagonal" degrees of free- 
dom of the metric tensor arise and the foregoing new 
approximate system reduces to the ordinary differen- 
tial equations (for each spatial point) for the variables 
J3 1 ,...,/3 d and fi d+l where j3 rf+1 = —(p. It is convenient 
to use the d + 1 -dimensional flat superspace with coordi- 
nates /3 1 , ...,/3 d+1 and correspondingly new indices A, B 
running over values l,...,d+l. The metric G A g m this 
superspace is 

G M dfdf = t (dP A ) 2 -ddl5 A ) 2 +(dl5 d+l ) 2 . 

A=l A=l V 7 

(20) 

The asymptotic dynamics for j3 -variables follows from 
the Lagrangian of the form similar to (14): 

Gi«rV = -i. 

(21) 

Again ([]) component of the Einstein equations gives 
additional condition to the equations of motion following 
from this Lagrangian: 

(22) 

Here j3 A = py 4 and time parameters T and % are de- 
fined by the evident generalization to the multidimen- 
sional spacetime of their definitions we used in case of 
4-dimensional gravity: dt = — y/detg a pdT, dx = p 2 dT . 
All functional parameters Cb(x a ) in general are positive. 
The cosmological singularity corresponds to the the limit 
p — > oo and in this limit potential term in Lagrangian can 
be considered as p -independent, asymptotically it vanish 
in the region of this space where Wb{j) > and is in- 
finite where w&(y) < O.The sum in the potential means 
summation over all relevant (dominating) impenetrable 
barriers located at hypersurfaces where w&(y) = in the 
hyperbolic <f -dimensional y-space. All w&(y) are linear 
functions on y therefore w&(y) = p~ l Wb(fi)- The free 
motion of j3 A (r) between the walls in the original d + 1 - 
dimensional /3 -superspace is projected onto a geodesic 
motion of y^{T) on hyperbolic J-dimensional y-space, 
i.e. to the motion between the corresponding projections 



of the original walls onto y-space. These geodesic mo- 
tions from time to time are interrupted by specular re- 
flections against the infinitely sharp hyperplanes w&(y) = 
0. These hyperplanes bound a region in y-space inside 
which a symbolic particle oscillates and the volume of 
this region, in spite of its non-compactness, is finite. The 
last property is of principle significance since it leads to 
the chaotic character of the oscillatory regime. 

Of course, one of the central point here is to find all the 
aforementioned dominant walls and corresponding "wall 
forms" Wfc(j3). This depends on the spacetime dimension 
and menu of p-forms. In papers [33, 7, 8] the detailed 
description of all possibilities for the all types of super- 
gravities (i.e., eleven-dimensional supergravity and those 
following from the known five types of the superstring 
models in ten-dimensional spacetime) can be found. It 
was shown that in all cases there is only 10 relevant walls 
governing the oscillatory dynamics. The large number of 
other walls need no consideration because they are lo- 
cated behind these principal ten and have no influence 
on the dynamics in the first approximation. The men- 
tioned above region in y-space where a particle oscil- 
late is called "billiard table" and collection of its bound- 
ing walls forms the so-called Coxeter crystallographic 
simplex, that is, in the cases under consideration, poly- 
hedron with 10 faces in 9-dimensional y-space with all 
dihedral angles between the faces equal to the numbers 
n/n where n belongs to some distinguished set of nat- 
ural numbers (or equal to infinity). This is very special 
geometric construction which (when combined with the 
specular laws of reflections against the faces) lead to the 
nontrivial huge symmetry hidden in the asymptotic struc- 
ture of spacetime near cosmological singularity which 
symmetry coexists, nevertheless, with chaoticity. 

DAMOUR-HENNAUX-NICOLAI 
HIDDEN SYMMETRY CONJECTURE 

The mathematical description of the symmetry we are 
talking about can be achieved in the following way. Con- 
sider the trajectories j3 A (t) of a particle moving between 
the walls Wb(j5) = in the original 10-dimensional j3- 
superspace with coordinates j3 A and metric G A g (20). 
These trajectories are null stright lines with respect to 
the Lorenzian metric G A g. Wall forms W&Q3) are lin- 
ear function on j3, that is Wb = w bA f5 A where the set of 
constants w hA depends on the choice of a supergravity 
model and on the type of the wall (index b) in the chosen 
model. We see that for each wall Wb = the constants 
w hA - represent components of the vector orthogonal to to 
this wall. We can imagine all these vectors (for differ- 
ent b) as arrows starting at the origin of the j3 -space. All 
these vectors have fixed finite norm G AB w b AWbB (G AB is 



inverse to G^g) and one can arrange the scalar products 
(w a • wi,) = G AB w a ^w b g for each supergravity model in 
the form of the matrix: 

A a b = 2 7 r (no summation in a). (23) 

The crucial point is that, independently of a supergrav- 
ity model, A a b is the Cartan matrix of indefinite type, i.e. 
with one negative principal value [8, 33, 34, 35]. Any 
Cartan matrix can be associated with some Lie algebra 
and particular matrix (23) corresponds to the so-called 
Lorenzian hyperbolic Kac-Moody algebra of the rank 
10. As was shown in [31] the particle's velocity v* 4 = 
d$ A jd% after the reflection from the wall w a ^fi A = 
changes according to the universal (i.e. again indepen- 
dent of the model) law: 

after _ = {^before - 2 ^^^ W i > ( 2 4) 

w% = G AB w a g, (no summation in a). 

This transformation is nothing else but the already men- 
tioned specular reflection of a particle by the wall or- 
thogonal to the vector w a ^. Now it is clear that one can 
formally identify the ten vectors w a ^ with the simple 
roots of the root system of Kac-Moody algebra, the walls 
w uaP A = with the Weyl hyperplanes orthogonal to the 
simple roots, the reflections (24) with the elements of 
the Weyl group of the root system and the region of j3- 
superspace bounded by the walls (where a particle oscil- 
lates) with the fundamental Weyl chamber. For the read- 
ers less familiar with all these notions of the theory of 
generalized Lie algebras (especially in application to the 
question under consideration) we can recommend the ex- 
haustive review [10] which is well written also from ped- 
agogical point of view. 

The manifestation of Lie algebra means that the cor- 
responding Lie symmetry group must somehow be hid- 
den in the system . The hidden symmetry conjecture 
[35, 36, 37] proposes that this symmetry might be inher- 
ent for the exact superstring theories (assuming that they 
exist) and not only for their classical low energy limits of 
their bosonic sectors in the vicinity to the cosmological 
singularity. The limiting structure near singularity should 
be considered just as an auxiliary instrument by means 
of which this symmetry is coming to light. As of now we 
have no comprehension where and how exactly the sym- 
metry would act (could be as a continuous infinite dimen- 
sional symmetry group of the exact Lagrangian permit- 
ting to transform the given solutions of the equations of 
motion to the new solutions). If true the hidden symmetry 
conjecture could create an impetus for the third revolu- 
tion in the development of the superstring theories. 
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